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Abstract 

In the presence of arbitrarily large magnetic fields, matter composed of electrons 
and nuclei was known to be unstable if a or Z is too large. Here we prove that matter 
is stable if a < 0.06 and Za 2 < 0.04. 

One of the remaining unsolved problems connected with the stability of matter is the 
inclusion of arbitrary magnetic fields. The model is a caricature of QED which invites 
speculations about stability of QED for large fine structure constant, a, but that is not our 
focus here and we refer to [0] for a discussion of these and related matters. The Hamiltonian 
for N electrons and K fixed nuclei of charge Ze with magnetic field B(x) = V x A(x), 
including the field energy, e J B 2 , is 

N 

H = Y J % + V c + e B(x) 2 d 3 x, (1) 
i=i J 

where T = [a ■ (p + A)} 2 — (p+ A) 2 + <r • B is the Pauli operator. The Coulomb energy is 

N K 

V c = -Z^2J2\xi - Rjl' 1 + \ x i~ x i\~ l 

i=l j=l l<i<j<N 

+z 2 Yl \Ri-Rj\~ 1 , (2) 
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with Rj being the coordinates of the nuclei and X{ the electron coordinates. The energy unit 
is 4 Rydbergs = 2mc 2 a 2 , a = e 2 /hc, length unit = h 2 /2me 2 and e = (Sna 2 )" 1 . Notice that 
a appears in ([[]) only through e. 

The negative particles, i.e., the electrons, are necessarily spin 1/2 fermions which, for 
generality, we assume to exist in q/2 flavors (e.g., q = 6). The ground state energy is denoted 
byE. 

Starting with the 1967 pioneering work of Dyson and Lenard we now understand stability 
for arbitrarily many electrons and nuclei, with B = 0, in the context of the nonrelativistic 
Schrodinger equation. Later it was extended to the "relativistic" Schrodinger equation in 
which p 2 /2m is replaced by (c 2 p 2 + m^c 4 ) 1 ^ 2 (see [0 for a review). These proofs also hold 
with the inclusion of a magnetic field coupled to the orbital motion of the electrons, i.e., 
p — > p + A, but no Zeeman cr • B term. 

Stability of matter has two meanings: (1) E is finite for arbitrary N and K; (2) E > 
—Ci(N + K) for some constant C\ independent of N, K and Rj. In the nonrelativistic case 
(2) holds. In the relativistic case, (1) implies (2) but (1) requires two conditions: Za < C 2 
and a < C3 with C*2 and C3 being universal constants, the best available values being in 
0, Theorems 1 and 2. The inclusion of B changes E, but the point is that while C\, C2, C3 
depend on q, they can be chosen to be independent of B. 

The situation changes dramatically when the magnetic moments of the electrons are 
allowed to interact with the magnetic field via the cr ■ B term, as in dH). The reason for this 
is simple: The Pauli operator T is nonnegative but it is much weaker than (p + A) 2 . Indeed, 
it can even have square integrable zero-modes 0, Tip = 0, for suitable A(x), which cause 
instability for large Za 2 . 

It is known 0] that without the field energy term, e J B 2 , in ([I]) arbitrarily large B fields 
can cause arbitrarily negative energies, E, even for hydrogen. The field energy, hopefully, 
stabilizes the situation, and our goal is to show that E is finite for (|l]), even after minimizing 
over all possible B fields and all possible Rj. 

One of our results on magnetic stability is: 

Theorem 1: The ground state energy of H satisfies 

E > -2.6 q 2 / 3 max{Q(Z) 2 , Q(5.7g) 2 }AT 1 / 3 ^ 2 / 3 , (3) 
with Q(t) =t + \plt + 2.2, provided that 

qZa 2 < 0.082 and qa < 0.12 . (4) 
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In (0) all the nuclear charges are set equal to Z . As far as stability is concerned this is 
no restriction |5| since the energy is concave in each charge Zj and hence stability holds in 
the "cube" {0 < Zj < Z}^ =1 if it holds when all Zj = Z. It also follows from this that E is a 
decreasing function of Z. Moreover, since e oc a~ 2 it follows that E is a decreasing function 
of a, a fact that will be important later. The form of (|3]) is the best possible for Z > 1, as 
we know from other studies ||. 

Our actual condition for stability given after fll8|) is rather complicated, but very much 
more general than (f|) — which is only representative. The results after fllq) show, e.g., that 
when a = 1/137, and q = 2, Z can be as large as 1050. The large values of Z and a are 
important because the comfortable distance of the critical values from the physical values 
Z < 92, a = 1/137 implies that the effect studied here is merely a small perturbation. 

Our proof of Theorem 1 will require a new technique — a running energy-scale renormal- 
ization of T. A by-product of this is a Lieb-Thirring type inequality for T: 

Theorem 2: If e\ < £2, ■ ■ ■ < are the negative eigenvalues of T — U, for a potential 

-U{x) < then 

EM < a iJ U{x) 5/2 d 3 x 

+ 6 7 (| B(x) 2 d 3 x y /4 (J f/(x) 4 d 3 xy /4 (5) 

for all < 7 < I, where a 7 = (2 3 / 2 /5)(l - 7)-% and b~ f = ^'^/^-^{l - 7)- 5/8 L 3 . 
We can take L 3 , defined below, to be 0.1156. 

More generally the second term in (JD can be replaced by (/ B^l 2 ) 1 li(f I/p) 1 ^, where 
p- 1 + q- 1 = 1. 

The investigation of this problem started in Jl| where type 2 stability was proved (for 
suitable Z,a and q = 2) for K = 1 and arbitrary iV (if (Z + l/4)a 12/7 < 0.15) or iV = 1 and 
arbitrary K (if Za 2 < 0.6 and a < 0.3) . The problem for general N and K was open for 9 
years and we present a surprisingly simple solution here. 

The bounds in (|J) on Za 2 and a are not artifacts. It is shown in and j|] that the 
zero-modes cause E = —00 when Za 2 > 11.11 for the "hydrogenic" atom, i.e., a single spin 
1/2 particle and one nucleus. If the number of nuclei is arbitrary, it is shown in [0] that there 
is collapse if a > 6.67, no matter how small Z is. Magnetic stability, like relativistic stability, 
implies a (Z-independent) bound on a. 

Prior to our work a proof of type 2 stability for (0) with Z = 1, q = 2, and some 
sufficiently small a was announced (unpublished) by C. Fefferman and sketched to one of us. 
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Our proof is unrelated to his, considerably simpler and, more importantly, gives physically 
realistic constants. 

We begin our analysis with the observation that length scaling considerations suggest 
that the key to understanding the stability problem is somehow to replace T, on each energy 
scale, e, by /xT/e, where \i is a fixed energy but e is variable. On energy scales e>/iwe can 
use the fact that T > to replace T by fiT/e without spoiling lower bounds. It might seem 
odd to replace T by something smaller, but what is really happening is that a ■ B is being 
partially controlled by [1 — /xe _1 ](p + A) 2 . The idea of replacing T by a fraction of T was 
also used in [JT], but no energy dependence was used there. 

We shall illustrate this concept by three calculations. The first, A, will establish magnetic 
stability by relating it to the stability of relativistic matter (see |3|, f| || ||). The second, B, 
will be the proof of Theorem 2. The third, C, will use essential parts of the second calculation 
and an electrostatic inequality proved in to prove magnetic stability without resorting to 
relativistic stability. 

A. Magnetic Stability from Relativistic Stability: We use stability of relativistic 
matter in the form proved in ||. From the corollary of Theorem 1 in with (3 = 0.5 we 
have, for any < qn < 0.032 and Zk < 1/vt, 

N 

ElPi + ^l+zcK^O. (6) 

i=i 

(Although Theorem 1 in 0] was stated only for it holds for \p + A\ because it relies 
only on the magnitude of the resolvent which only gets smaller when A is not zero. That is 
, ||p + A\~ s (x,y)\ < \\p\~ s (x,y)\ for each s > and x,y in i? 3 . This follows at once from 
a similar bound on the heat kernel {exp[— t(p + A) 2 ]}(x,y) which, in turn, follows from its 
representation as a path integral. This was pointed out in [|, |i~0fl . Only the resolvent powers 
\p + A\~ s enter the proof of Theorem 1 in 0.) 

Using dB|), if is bounded below by H — Y$Li hi where h is the one-body operator h = 
T — K^lp + A\. Thus, E is bounded below by e J B 2 + E N} where E N = qY^jHi £ j an d 
£\ < ^2 ^ ••• are the eigenvalues of h. For e > 0, let N_ e (h) be the number of eigenvalues of 
h less than or equal to — e. Choose fi > and note that 

E N > —N/jl -q N_ e (h)de . (7) 
Jn 

The crucial step in our proof is noting that the positivity of the operator T implies that 
T > fiT/e when e > /x. Thus, T > /ie _1 T > /ie _1 (p + A) 2 — \ie~ x B{x) when e > ji. 
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By Schwarz's inequality, k 1 \p + A\ < (l/3)e l K 2 {p + A) 2 + (3e/4) and hence if we set 
ji = (4/3)k -2 we obtain 

h > e- X K- 2 (p + A) 2 - (A/3)e- 1 n- 2 B(x) - (3e/4) = /i e . 

Thus, N_ e (h) < N^ e (h e ) and this can be estimated by the Cwikel-Lieb-Rozenblum (CLR) 
bound 0, i.e., N_ e ((p + A) 2 - U(x)) < L 3 J[U(x) - e] s / 2 d 3 x where [a]+ = max(a, 0) and 
L 3 = 0.1156. In our case: 

3/2 



N. e {h e ) <L 3 J 



AB(x) e 2 K 2 



d s x . (8) 
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Inserting this bound in (0), a simple calculation yields 

E N > -Nfi - (2n/3)qK- 1 L 3 J B(x) 2 d 3 x . 

We choose k so that the field energy terms are non-negative, i.e., k > (16ir 2 /3)L 3 a 2 q = 
6.1a 2 q. We conclude, by (|6|), that magnetic stability holds if 

qa < 0.071 and qZa 2 < 0.052 . (9) 

For q = 2, the first condition is a < 1/28. For q = 2 and a = 1/137 stability occurs if 
Z < 490. 

Assuming (^) holds, we then use (^|) and choose k = min{0.0315g _1 , (irZ)~ 1 }. Our lower 
bound on the ground state energy per electron, by this method, is then — /j, = — (4/3)/t~ 2 = 
-max{1345g 2 , 13.2Z 2 }. 

Remark: We used the CLR bound in (||). Since the derivation of this bound is not 
elementary the reader might wish to use an easier to derive bound — at the cost of worsening 
the final constants. A useful substitute is 



N_ e < 0.1054e- 1/4 y" [U(x) - e/2} 7 / 4 d 3 x 



which is in (2.8) of [|I^] and which can be derived by means originally employed for the 
Lieb-Thirring inequality. This same remark also applies to our other calculations below. 



B. The Lieb-Thirring Inequality: As before we note that Yl £ i — ~ Jo° N_ e (T — U)de. 
We write J °° = Jq + J^°. The parameter \i will be optimized below. Noting that T > 
(p + A) 2 — B(x) and applying the CLR bound in the same fashion as before to Jq 1 yields 

L 3 J \B{x) + U(x) - e]+ /2 d 3 xde . (10) 
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In we replace T by the lower bound fie 1 [(p + A) 2 — B(x)} and obtain AL e (T — U)< 
./V_ e (/ie -1 [(p + A) 2 — 5] — U). A further application of the CLR inequality yields the bound 

TOO 

on /„ 

/•oo r _ ,_ 

(11) 



L 3 j j [B(x) + (e/fi)U(x) - e 2 //i]+ /2 dxde . 
It is easy to see that for any < 7 < 1 the integrand in (|T0|) is bounded above by 

V2 {[B(x) - 7 e 2 /^]f + P{x) - (1 - 7 )e]f ) . 
Treating the integrand in (|ll]) in a similar fashion and combining the inequalities we find 

E N < ^3 / (jTW) - 7 e 2 /^]f de 
+ (1 - 7 ) e ]fde 

{e/n)U{x) - (1 -7)e 2 //il 3/2 de)d 3 x . 



+ 

After extending the last two integrals to / °°, a straightforward computation yields 

2 3V/ 2 



EN< ^/{^ 



7) 



U(x 



,5/2 



+ 



I67V2 



fix 



3tt 



f^- 3/2 (l-7)- 5/2 ^) 4 Kx. 



Optimizing over /i yields (|5|). 

To prove the more general form of (|5|) replace /ze -1 by (/ie^ 1 ) 5 , where s = 2p/3 — 5/3. 



C. Proof of Theorem 1: We turn now to our third illustration of the concept of running 
energy scale and prove the stability directly, not relating it to the relativistic problem. By this 
method we get the correct dependence of the ground state energy on Z and also somewhat 
better critical constants than in (||). 

Following P| we first replace the Coulomb potential by a single particle potential in (|T^) 
below. We break up -R 3 into Voronoi cells defined by the nuclear locations, i.e., Tj = {x : 
\x — Rj\ < \x — Rk\ for all k} is the j-th Voronoi cell. Each Tj contains a ball centered at Rj 
with radius Dj = mm{\Rj — R k \ : j 7^ k}/2. 

The following bound on V c is proved in ||: Choose some < A < 1. Then 

N K y2 

K>-£w^) + E-— , (12) 

i=i j=i 
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where W(x) = Z\x — Rj\ 1 + Fj(x) for x G Tj with Fj(x) defined by 

{2D j )~ 1 {l - Dj 2 \x - Rj\ 2 Y l for \x -Rj\< XDj 



2Z+l/2)\x-R j \- 1 for \x-R j \>XD j . 

The point about this inequality is that the potential W has the same singularity near each 
nucleus as V c has, and that the rightmost term in fljjD is repulsive. This term will be 
responsible for stabilizing the system. 

The problem is thus reduced to obtaining a lower bound on qJ2' e j, where J2' £ j is th e 
sum of the first [N/q] negative eigenvalues of of T — W. Note that Theorem 2 cannot be 
applied directly to this problem, since W is neither integrable to the power 5/2 nor to the 
power 4. Instead we have to do the calculations directly. 

For v > (a number that is chosen later) set W v (x) = (W(x) — u) + and note that 
W{x) -v< W v {x). Then, as in £j > -Nv - g/ °°A^_ e (T - W u )de. Again, 

/ N_ e (T-W v )de < / N- e (T — W u )de 
Jo Jo 

roo 

+ / N_ e {iie~ x 7 - W)de , (13) 
J n 

where we have replaced W v (x) by W(x) in the second term. Applying the CLR bound to 
the first expression on the right side we obtain L3 / Jq[B(x) + W u (x) — e]!^ 2 ded 3 a; which can 
be bounded, as in part B, by 



M - - • 7e 2 i3/2 



V2L 3 I\ J^[B(x) - de 



+hi- 1 y 1 w l/ (xf 2 \d 3 x , (14) 
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for any < 7 < 1. 

The difficulty in dominating the second term in ([13]) comes from the Coulomb singularity 
of W(x) which is not fourth power integrable. The singularity can be controlled using the fol- 
lowing operator inequality which follows from the diamagnetic inequality / \ {p + A)ip\ 2 d 3 x > 
J|p|^|| 2 d 3 x and Lemma 2a on p. 708 of [[L3|]. 

- 1 if |x| < R 
if Id > R 



(P+ A?-zi\4>-\ z ;iit (mzR 

\Zj\X\ , 



Choose R = XDj and write (p + A) 2 = (3(p + A) 2 + (1 - (3)(p + A) 2 for some < f3 < 1. 
Then, by scaling, 

(/i/e)T -W v > (/i/e)(l -/3)(p + A) 2 - ^/e)B - W , 
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where W(x, e) = Gj(x, e) + Fj(x) for x G Tj with Gj(x, e) defined by 

(Z 2 e/Apfi) + 3Z/(2AD j ) for \x - Rj\ < XDj 
Z\x — Rj\~ for |x — Rj\ > XDj . 

Note that W depends on e. 

Again, as in part B, we can use the CLR bound on the second term in ([13|) to obtain 
[when 1 - 7 > Z 2 /(A/3fi)\ 

V2 L 3 (l - /3)- 3/2 / { j™\B{x) - 7 e 2 //i]fde 

+/^ 3/2 J™[eW(x, e) - (1 - 7 )e 2 ]+ /2 de|d 3 x . (15) 

First we compute the last integral in fllS]) , which is 

-K" /• POO 

V / / [eG,-(x, e) + eF,(x) - (1 - 7 )e 2 ] 3 / 2 ded 3 £ . 
i=1 J \, Jo 

Now split the Tj integral into an inner integral \x— Rj\ < XDj and an outer integral \x — Rj\ > 
XDj. The inner integral yields, using the definitions of Gj and Fj, 



32 v ' 4/V 



2(1 -r 2 ) 

r 2 drD/ 1 . (16) 

To bound the outer integral from above we replace Tj by -R 3 and get 



3Z 
+ 2A 



(3ir/32)(l-7r 5/2 (VZ+ V 'l/2) 8 (A, Dj )- 1 . (17) 

Combining (p^)-(|17|) we find that the sum of the negative eigenvalues of T — W v is bounded 
below by 

K 

-a J W„{x) 5/2 d 3 x - b J B(x) 2 d 3 x - cY^Dj 1 . (IN) 
Here a = q{2V2/b)L 3 {l - 7)- 1 , 

b = ^i 3 (l-^)- 3/2 W7) 1/2 
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3=1 



9 32 31 P) ^ I A(l- 7 ) 5 / 2 

V2 „x 4 

r 2 dr [. 

To simplify the stability condition we have artificially increased the bounds by recalling that 
q > 2 and twice replacing 1/2 by q/4 in the definition of c. We choose f3 = 1/8, 7 = 1/2, 
A = 8/9 and /i so that b = (871a 2 )" 1 . The stability condition c < Z 2 /8 [see (|T2"D] now depends 
on/?/ on the 2 parameters X = qZa 2 and F = qa. A straightforward, but lengthy calculation 
shows that the stability condition holds if X = X Q = 0.082 and Y = Y = 0.12. The condition 
is monotone in Y, so it holds for X = Xo, Y < Y$. Although our condition does not hold for 
all X < Xq, Y < Yq we can use the Z-monotonicity of E to conclude stability in this range; 
this proves (|]). With the same values of 7 and A and with q = 2 the values Z = 1050, 
a = 1/137 also give stability. 

To derive (|^), note that W(x) < Q\x— -Rj| _1 for x G Tj. Using this bound and replacing Tj 
by i? 3 , one easily obtains —^2-n 2 L^qKQ^v^l 2 — X v as a lower bound on the —a J W^ 2 term 
in (P~5|) . Optimizing over v yields (^) when X = X , Y < Y . In this case, Z > Z = 5.7q. 
If X < X , Y <Y and Z > Z we get a lower bound on E by increasing a until X = X Q , 
y < this yields (||) with Q = Q(Z). Otherwise, with Z < Z , we use the Z-monotonicity 
of E to conclude (|3|) with Q = Q(5.7q). 

This work was partially supported by NSF grants PHY90-19433-A04 (E.H.L.), DMS92- 
07703 (M.L.) and DMS92-03829 (J.P.S.). 



+ (1-7- 



AR,,' 



-5/2 



./in 



7r{ + oT. 



References 

[1] J. Frohlich, E. Lieb and M. Loss, Commun. Math. Phys. 104, 251 (1986). 
[2] E. Lieb, Bull. Amer. Math. Soc. 22, 1 (1990). 

[3] E. Lieb and H.-T. Yau, Commun. Math. Phys. 118, 177 (1988); Phys. Rev. Lett. 61, 
1695 (1988). 

[4] M. Loss and H.-T. Yau, Commun. Math. Phys. 104, 283 (1986). 

[5] J. Avron, I. Herbst and B. Simon, Duke Math. J. 45, 847 (1978); Com- 
mun. Math. Phys. 79, 529 (1981). 



9 



[6] I. Daubechies and E. Lieb, Commun. Math. Phys. 90, 497 (1983). 

[7] E. Lieb and M. Loss, Commun. Math. Phys. 104, 271 (1986). 

[8] J. Conlon, Commun. Math. Phys. 94, 439 (1984). 

[9] C. Fefferman and R. de la Llave, Rev. Math. Iberoamericana 2,119 (1986). 

[10] J. Combes, R. Schrader and R. Seiler, Ann. Phys. (NY) 111, 1 (1978). 

[11] E. Lieb, Proc. Amer. Math. Soc. Symposia in Pure Math. 36, 241 (1980). 

[12] E. Lieb and W. Thirring in Studies in Mathematical Physics, edited by E. H. Lieb, 
B. Simon, A. Wightman (Princeton Univ. Press, Princeton 1976), p. 269. 

[13] A. Lenard and F. Dyson, J. Math. Phys. 9, 698 (1968). 



10 



